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ABSTRACT 

We assess the effects of simulated active galactic nuclei (AGNs) on the colour and 
morphology measurements of their host galaxies. To test the morphology measure- 
\ ments, we select a sample of galaxies not known to host AGNs and add a series of 

point sources scaled to represent specified fractions of the observed V band light de- 
tected from the resulting systems; we then compare morphology measurements of the 
simulated systems to measurements of the original galaxies. AGN contributions >20 
£Nl ■ per cent bias most of the morphology measurements tested, though the extent of the 

<0 apparent bias depends on the morphological characteristics of the original galaxies. 

We test colour measurements by adding to non-AGN galaxy spectra a quasar spec- 
trum scaled to contribute specified fractions of the rest-frame B band light detected 
from the resulting systems. A quasar fraction of 5 per cent can move the NUV— r 
colour of an elliptical galaxy from the UV-optical red sequence to the green valley, 
and 20 per cent can move it into the blue cloud. Combining the colour and morphology 
results, we find that a galaxy/AGN system with an AGN contribution >20 per cent 
may appear bluer and more bulge-dominated than the underlying galaxy. We con- 
clude that (1) bulge-dominated, E/SO/Sa, and early-type morphology classifications 
are accurate for red AGN host galaxies and may be accurate for blue host galaxies, 
unless the AGN manifests itself as a well-defined point source; and (2) although highly 
unobscured AGNs, such as the quasar used for our experiments, can significantly bias 
the measured colours of AGN host galaxies, it is possible to identify such systems 
by examining optical images of the hosts for the presence of a point source and/or 
measuring the level of nuclear obscuration. 

Key words: galaxies: active - galaxies: fundamental parameters - galaxies: nuclei - 
galaxies: structure. 



1 INTRODUCTION 

A galaxy's colours and morphological characteristics supply 
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example, optical colours allow one to estimate the ages of 
observed stellar populations (e.g., Kennicutt 1998), while 
UV-optical colours are remarkably sensitive to recent star 
formation (e.g., Wyder et al. 2007). Morphological char- 
acteristics provide information about a galaxy's interaction 
history, and when combined with the recent star formation 
history they can be used to estimate the amount of time that 
has passed since the most recent galaxy merger (e.g., Lotz 
et al. 2008b) . This is particularly useful for investigating po- 
tential causes of significant black hole growth (active galac- 
tic nuclei; AGNs), and possible connections between galaxy 
interactions and AGNs have been explored using both the- 
oretical and observational techniques. 

Many current simulations and semi-analytic models of 
galaxy interactions and evolution incorporate the effects 
of central supermassive black holes in the initial galax- 
ies (e.g., Kauffmann & Haehnelt 2000; Hopkins et al. 
2005a,b, 2008a,b; Croton et al. 2006; Ciotti & Ostriker 2007; 
Somerville et al. 2008). The results of these simulations in- 
variably predict that the black holes become active dur- 
ing some stage of a merger event between gas-rich galax- 
ies of similar masses (a 'major merger', typically defined as 
a merger between galaxies having a mass ratio between 1:1 
and 3:1). Most of these authors additionally found that some 
form of radiative output (or 'feedback') from the AGN sig- 
nificantly decreased the star formation rate in the galaxy. 
Hopkins et al. (2005a,b, 2008a,b) described a form of AGN 
feedback that forcefully removes gas from the nuclear re- 
gions of a galaxy, after significant star formation and ac- 
cretion onto the black hole has used up much of the gas. 
Croton et al. (2006) presented a complementary scenario, a 
'radio mode', in which AGN feedback heats gas in the nu- 
clear regions to temperatures sufficient to stop star forma- 
tion, without necessarily removing the gas from the galaxy. 

Schawinski et al. (2007) and Georgakakis et al. (2008) 
described recent observational support for the scenarios de- 
rived from simulations (such as those summarized above) 
and semi-analytic models (e.g, Somerville et al. 2008). After 
visually examining ~16,000 low-redshift early-type galaxies 
and spectroscopically classifying each as experiencing star 
formation, black hole growth, a combination of the two, or 
neither, Schawinski et al. (2007) presented evidence that ra- 
diative output from an AGN may cause its host galaxy to 
transition from a star-forming system to a quiescent system. 
They reported that the optical colours of the emission-line 
systems (20 per cent of their sample) are typically bluer 
than the quiescent systems, and the magnitude of the off- 
set from the red sequence toward the blue cloud depends 
on the type of activity observed. Most of the systems domi- 
nated by emission from low-ionization nuclear emission line 
regions (LINERs) joined the inactive galaxies populating the 
red sequence, while the remaining LINER and Seyfert host 
galaxies typically exhibited colours on the red edge of the 
green valley. Early-type systems undergoing star formation 
were found in the blue cloud, and the systems exhibiting 
emission lines characteristic of both an AGN and star for- 
mation typically populated the green valley. 

Georgakakis et al. (2008) stacked the X-ray emissions 
of three samples of galaxies at z ~ 0.7. The samples were 
comprised of galaxies with optical colours placing them on 
the red sequence, the green valley, or the blue cloud, and 
were limited to systems from which X-rays had not been de- 



tected at the limits of the survey. After stacking, they found 
significant X-ray detections from the red sequence and green 
valley samples, but not the blue cloud, indicating the pres- 
ence of a significant population of obscured AGNs among 
red sequence and green valley galaxies. They concluded that 
accretion onto the central supermassive black hole may con- 
tinue after the star formation has ceased, and thus if AGN 
feedback is responsible for the cessation of star formation in 
these galaxies, then the method by which this occurs does 
not completely remove gas from the nuclear regions. 

In order to accurately compare theoretical predictions 
to conclusions based on observational studies, it is essential 
that we understand how the presence of an AGN may bias 
the colour and morphology measurements of its host galaxy. 
Morphologically, an unobscured AGN can be approximated 
as a point or as a highly concentrated object. Thus morphol- 
ogy measurements would be expected to experience a bias 
toward more bulge-like or spheroidal systems; the signifi- 
cance of this bias should depend strongly upon the balance 
between AGN and galactic light contributions. If the bias is 
significant and common, then we would need to reconsider 
the current classification of most AGN host galaxies as el- 
liptical or early-type (e.g., Grogin et al. 2003, 2005; Sanchez 
et al. 2004; Pierce et al. 2007), with occasional morphologi- 
cal evidence of recent interactions (e.g., Bennert et al. 2008; 
Urrutia, Lacy & Becker 2008). 

One plausible explanation for the elliptical classifica- 
tion is that AGNs tend to be observed in galaxies that have 
recently undergone a merger. However, if the morphology 
measurements are significantly biased, then many of the ap- 
parently post-merger systems hosting AGNs may in fact be 
disc-dominated, necessitating an alternate explanation for 
the black hole growth. A potential bias toward elliptical mor- 
phologies may also cause an underestimate of the merger 
rate, as determined by morphology measurements such as 
G-M20 (Lotz et al. 2008a). Thus, the proper classification of 
AGN host galaxy morphologies offers important clues as to 
when and how significant black hole growth is activated. 

In a recent study of AGN host galaxy morphologies, 
Gabor et al. (2009) accounted for the presence of a point 
source component in the following two ways: (1) using GAL- 
FIT (Peng et al. 2002) to measure two-dimensional surface 
brightness profiles, Gabor et al. (2009) utilized the option 
to include a point spread function (PSF) as a model for the 
nuclear point source in the fitting procedure; (2) they also 
measured the concentration and asymmetry of the galaxies 
both before and after subtracting a modeled point source 
that had been fitted to the AGN. They found that both 
concentration and asymmetry measurements are affected by 
the presence (or subtraction) of a point source. 

For the current study, we do not attempt to model a 
point source in the galaxies analyzed, but instead focus on 
understanding the effect that the presence of a point source 
may have on common morphology measurements. In AGN 
host galaxies, it is often unclear what fraction of the light 
emanates from the AGN and how much originates in star 
formation regions. A model of the nuclear light may thus 
overestimate the AGN contribution, and morphology mea- 
surements (such as G and M20) of the residual from a sub- 
sequent subtraction may experience a different set of biases 
than those caused by the actual point source. An investiga- 
tion of such effects is beyond the scope of this paper. 
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Kim et al. (2008) emphasized the importance of using 
an appropriate PSF for studies similar to that described 
by Gabor et al. (2009), and they found that PSF charac- 
teristics can vary strongly with time, location on a camera 
detector, and level of pixel sampling. Mismatches between 
the PSF appropriate for a particular AGN and the PSF(s) 
determined for a set of observations may cause inaccurate 
measurements of the host galaxy and/or AGN fluxes. How- 
ever, in § 13.21 we find that the choice of PSF is negligible for 
the present study. 

Measured colours may also be biased by light con- 
tributed by an AGN. Because optical and UV light from 
an AGN is typically much bluer than the light from its host 
galaxy (e.g., Kinney et al. 1996), a significant contribution 
from the AGN to the observed light may cause an intrinsi- 
cally red AGN host galaxy to appear blue or 'green', compli- 
cating or invalidating the use of colour to describe the host 
galaxy's star formation history. Sanchez et al. (2004) mea- 
sured the colours and morphologies of 15 AGN host galaxies 
at z ~ 1, after first subtracting a PSF individually deter- 
mined for each AGN. Upon comparing the resulting nuclear 
colours to the host galaxy colours, they found no adverse 
affects from the PSF subtraction. Though a few of their 
AGN host galaxies exhibit red colours, most show colours 
that are bluer than typical early-type galaxies, indicative 
of young stellar populations and consistent with the results 
from Schawinski et al. (2007) described above. 

The current paper investigates the biases caused by an 
AGN on five morphology measurements and two colour mea- 
surements. Four of the five morphology measurements con- 
sidered in the current work were initially calibrated without 
accounting for possible effects due to the presence of active 
nuclei (Lotz, Primack & Madau 2004). The fifth measure- 
ment technique includes an option to account for the pres- 
ence of a point source (Peng et al. 2002). AGN contamina- 
tion of UV-optical colours was recently tested by Kauffmann 
et al. (2007), who reported that the amount of AGN light po- 
tentially contaminating measured UV-optical colours should 
be minimal for low-redshift (z < 0.07), optical spectra- 
selected AGNs. Seeking to more fully understand the effect 
of higher-luminosity and higher-redshift (z ~ 0.4) AGNs on 
morphology and colour measurements, this paper presents 
an examination of the extent to which an AGN similar to a 
quasar or a Type-1 Seyfert galaxy may affect measurements 
of the Gini coefficient G, the second order moment of the 
brightest 20 per cent of the galaxy's flux M20, concentration 
C, asymmetry A, the Sersic index n, and optical (17 — B) 
and UV-optical (NUV— r) colours. 

We begin with a description of the two sets of data used 
for our experiments (§[21)5 followed by descriptions of the ef- 
fects of an AGN on the measured morphology (§ [3| and 
colour (§ [4]) of its host galaxy. The two sets of results are 
combined in § [SI where we present a potential track of an 
unobscured AGN in colour-morphology space, as a function 
of the fraction of B band light contributed by the AGN. Fi- 
nally, in §[U]we discuss and summarize our results. Through- 
out, {h, f^A, £Im} = {0.7, 0.7, 0.3}, and AB magnitudes are 
used, unless otherwise noted. 



2 DATA 

2.1 AGN host galaxy images 

The present study uses data available from the All- 
wavelength Extended Groth Strip International Survey 
(AEGIS; Davis et al. 2007, and references therein). In par- 
ticular, we use optical images taken with the Hubble Space 
Telescope/ Advanced Camera for Surveys (HST/ ACS), op- 
tical spectra taken with the DEIMOS spectrograph (Faber 
et al. 2003) as part of the DEEP2 Redshift Survey (Davis 
et al. 2003, 2007) and with the MMT Observatory (Coil et 
al. 2009), deep Chandra X-ray Observatory images (Laird 
et al. 2009), and Very Large Array (VLA) images taken at 
1.4 GHz (Ivison et al. 2007). The galaxies used for our first 
experiment are drawn from a region of AEGIS for which all 
of these data are available. 

We select galaxies that have spectroscopic redshifts 
0.2 < z < 0.6, so that all the analyses can be performed 
using HST/ ACS V band images, which approximately cor- 
respond to the rest-frame B band at these redshifts. The 
HST/ ACS V band images analysed here were taken with the 
F606W filter, using single-orbit exposures and four dithers. 
The resulting images have a PSF of ~ 0.1" and a scale of 
0.03"/pixel; see Lotz et al. (2008a) for details regarding the 
reduction of the HST/ ACS images. 

Galaxies that have X-ray luminosities L2-10 kcv > 10 42 
erg s" 1 , radio powers Pi. 4 GHz > I0 1 ' 2z+23 2 W Hz" 1 , and/or 
have been spectrally identified as Seyfert galaxies, LINERs, 
or broad-line AGNs are excluded from the sample in order 
to control the effect of an AGN0 We intentionally exclude 
AGNs from our sample because the current study aims to 
determine the effect that an AGN might have on colour and 
morphology measurements of its host galaxy, and including 
AGNs would severely complicate these efforts. 

A sample of galaxies meeting these criteria is created 
from 29 of the 63 HST/ ACS tiles, each containing between 
12 and 34 eligible galaxies. This sample of 622 galaxies is 
later refined by imposing criteria designed to assure the re- 
liability of the individual morphology measurements (§ 13.10 . 
Because each of the individual and combined morphology 
measurements require different criteria, the final samples 
used for each morphology analysis differ slightly; additional 
descriptions of each sample are provided with the results. 



2.2 AGN spectral template 

Spectral templates are available for a variety of AGN types, 
including quasars (e.g., Francis et al. 1991), LINERs, and 
Seyfert galaxies (e.g., de Bruyn & Sargent 1978). Compared 
to a quasar, lower luminosity and/or more heavily obscured 
AGNs are expected to less significantly bias the measured 
host galaxy colours, although some contamination will be 
unavoidable for all but the most heavily obscured AGNs. 
One goal of the current work is to determine the level of this 
bias in AGN host galaxies. Ideally, we would determine this 
for all AGN types, however, the intervening galactic matter 



1 Pierce (2009) described in greater detail the identification of X- 
ray, radio, and optical spectra-selected AGNs in the AEGIS and 
GOODS-N regions. 
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Figure 1. Non-AGN (solid blue curves) and quasar (dot-dashed 
red curves) spectral templates. Gray shaded regions represent the 
effective widths of the four passbands indicated. Note that the U 
and B passbands do not overlap; the solid vertical line indicates 
where they split at A = 3980A. 



Observatory (CTIO), and the Kitt Peak National Observa- 
tory. Optical spectra observed at the CTIO were measured 
with slits 10" x 20", designed to match the area of the IUE 
aperture. Table 1 of Kinney et al. (1996) lists properties of 
the individual galaxies observed, and their Table 2 specifies 
which galaxies are used for each template. 

To create the spectral templates, Kinney et al. (1996) 
shifted each spectrum to its rest frame, corrected for Galac- 
tic extinction, and then averaged the UV and optical spec- 
tra within each galaxy morphology group (e.g., Sb). Optical 
spectra were averaged separately from UV spectra, but sim- 
ilar methods were employed for both averaging processes, 
and then the two averaged spectra were combined to pro- 
duce the final templates (shown in Figure 2 of their paper); 
the templates used for the current paper are shown in Fig- 
ure [T] Kinney et al. (1996) compared their non-AGN galaxy 
templates to results from Kennicutt (1992) and noted that 
their Sb template does not sufficiently represent young stel- 
lar populations in the disc. This suggests that colours de- 
rived from the Sb galaxy template may be too red. However, 
this should not significantly affect the work reported here. 



present in obscured AGNs and the significant light contribu- 
tions from the host galaxies of lower luminosity AGNs will 
contaminate the observed AGN light, biasing the measured 
colour of the AGN itself. Since this is the reverse of the 
problem that we want to address with the current study, 
we select a quasar template in order to avoid the poten- 
tially significant light contamination from the host galaxy. 
Although quasars are among the less common and more ex- 
treme AGN types at redshifts z ~ 1, the current work pro- 
vides upper limits to the biases in the measured AGN host 
galaxy colours. 

Francis et al. (1991) prepared a composite quasar spec- 
tral template based on more than 700 individual quasar and 
AGN spectra from the Large Bright Quasar Survey. The 
observations covered at least 3200-6800A at a resolution of 
approximately 6A and a S/N= 12. After correcting for red- 
shifts, between 200 and 500 AGNs contribute to the compos- 
ite spectrum at rest-frame wavelengths Ao = 1200-4200A, 
while up to 200 contribute at A = 800- 1200 A and A = 
4200-6000A. 

Figure Q] shows the quasar composite spectrum along 
with the spectra of the three non-AGN galaxy templates 
described below; note that although the quasar template 
does not extend all the way across the r band, we do not 
expect this to be problematic, due to the low contribution 
expected from the quasar in the r band, relative to the U 
and NUV bands. The most prominent feature of the quasar 
spectral template, as shown in Figure [1] is the broad Mgll 
line at A = 2798A. H/3 (A = 486lA) and [OIII] (AA = 4959, 
5007A) lines are also apparent. [Figure 2 of Francis et al. 
(1991) shows the quasar spectrum in greater detail.] 

2.3 Non-AGN galaxy spectral templates 

Kinney et al. (1996) presented spectral templates of sev- 
eral galaxy types - including elliptical, Sb, and starburst 
galaxies - based on observations using the International Ul- 
traviolet Explorer (IUE), the Cerro Tololo Inter- American 



3 EXPERIMENT #1: THE EFFECT OF AN 
AGN ON MORPHOLOGY MEASUREMENTS 

3.1 Classifying galaxy morphologies 

Galaxy morphologies depend upon the rest-frame wave- 
lengths of the images analysed. At any given redshift, the 
observed wavelengths that correspond to the rest-frame B 
band are highly sensitive to recent star formation and young 
stellar populations; this distinction is particularly relevant 
for detection of recent or on-going galaxy interactions. Red- 
der bands correspond to older stellar populations and may 
miss indications of recent activity, while bluer bands (e.g., 
NUV) will miss stellar populations that are no longer dom- 
inated by the youngest stars (e.g., O and B). At the red- 
shifts used for the analyses described here (0.2 < z < 0.6), 
rest-frame B band light approximately corresponds to the 
observed HST/ ACS V band images. 

Code written by J. Lotz (e.g., Lotz et al. 2004) simul- 
taneously measures A, C, G, and M20, along with several 
additional quantities, from optical galaxy images. Lotz et 
al. (2004, 2008a) described a morphological classification 
system, using G and M20, that is particularly sensitive to 
interacting galaxies and is reliable to redshifts z ~ 1. In ad- 
dition, Lotz et al. (2008b) have shown that A and G may 
be more effective than C and A at separating interacting 
galaxies from non-interacting galaxies. 

Concentration and asymmetry (Conselice 2003) have 
been widely used by many authors to describe AGN host 
galaxies (e.g., Grogin et al. 2003, 2005); Pierce et al. (2007) 
instead used G and M20 for the same task. The light pro- 
files of AEGIS, GOODS, GEMS, and COSMOS galaxies (~ 
490,000 sources) have been measured by Griffith et al. (in 
preparation), using the GALAPAGOS code written by M. 
Barden (Barden et al. 2010), which implements the GALFIT 
program (Peng et al. 2002) to determine Sersic indices. 



2 http : //astro-staff. uibk. ac. at / ~m. barden / galapagos / 
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3.1.1 Asymmetry 

Measuring the asymmetry A of a galaxy is a three-step pro- 
cess, beginning with a rotation of the galaxy image (which 
includes all galaxy pixels within 1.5 r e n, where r e n is the el- 
liptical Petrosian radius; Petrosian 1976) by 180° about its 
central axis; in practice, the central axis is determined by 
minimizing A (Lotz et al. 2004). Next the rotated image is 
subtracted from the original image, and finally, the sum of 
the absolute value of the residual flux is compared to the 
flux of the original galaxy (Conselice et al. 2000; Conselice 
2003). This process is summarized by the following formula 
(Conselice et al. 2008): 



.4 



S|Jp — Ji8oj 

Sl/ol 



E|g - B 1 
ELBo] 



(1) 



where the subscripts refer to the rotation of the image, and 
I and B refer to the original image and the background, re- 
spectively. Highly symmetrical galaxies, such as undisturbed 
ellipticals or discs, leave minimal residuals, and A is small; 
the significant residuals of disturbed and irregular galaxies 
correspond to higher values of A. Galaxies with < A < 1 
are typical. 

3.1.2 Concentration 

Concentration C compares the radius enclosing a large frac- 
tion of a galaxy's flux to the radius enclosing a smaller frac- 
tion of its flux. For this work, we used the following definition 
(Conselice 2003): 



C = 5 x log(r 80 /r- 2 o), 



(2) 



where r n represents the radius enclosing n per cent of the 
flux within 1.5 r c n of the galaxy's centre. High values of C 
correspond to highly concentrated systems, such as ellipti- 
cal galaxies, and the value of C decreases for more diffuse 
galaxies, generally including disc galaxies and many inter- 
acting systems. Concentrations 2 < C < 5 are typical values 
for observed galaxies. 

3.1.3 The Gini coefficient 

Abraham et al. (2003) were the first to use G to describe 
the distribution of light among the pixels associated with 
a galaxy, though it was initially introduced as an economic 
statistic used to characterize the distribution of wealth in 
a society (Gini 1912). Lotz et al. (2004, 2008a) describe in 
detail the measurement of G and M20; here we provide brief 
summaries. 

If n is the number of pixels assigned to a galaxy and 
fi represents the flux from pixel i (ordered such that fi 
increases with the pixel index), then we can use the methods 
described by Lotz et al. (2004, 2008a) and measure G using 



G 



\f\n(n - 1) 



1)1/* I 



(3) 



High values of G correspond to galaxies in which the light is 
distributed among a small fraction of the associated pixels, 
such as a galaxy featuring a single bright nucleus or a galaxy 
that features multiple bright nuclei; this latter characteristic 
contributes to the sensitivity to interacting galaxies associ- 
ated with the G-M20 and G-A classifications. Low values of 



G correspond to a homogeneous distribution of light among 
the galaxy's pixels, such as would be found for a smooth disc 
galaxy. Gini measurements of 0.3 < G < 0.7 are common 
for observed galaxies. 



3.1.4 M 20 

The value of M20 depends on the spatial distribution of the 
light in a galaxy, relative to the centre of the galaxy. Lotz 
et al. (2004) defined M 20 as follows: 



M 2 . 



lo Sio 
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where 
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Mtot = £ M i = f^ Xi ~ + ^ ~~ 



Ve 



(4) 



(5) 



Here, the fluxes fi are arranged in order of decreasing flux, 
and (x c , y c ) is the galaxy centre, as determined by mini- 
mizing Mtot • The presence of multiple bright nuclei that are 
spatially separated from a galaxy's centre increases the value 
of M20 toward —1.0. However, if the brightest 20 per cent 
of the galaxy's light is positioned very close to the galactic 
centre, as would be the case for an undisturbed elliptical 
galaxy, the value of M20 approaches —2.5 or —3.0. Typical 
observed galaxies have — 0.5 > A/20 > —2.5. 

3.1.5 Reliability criteria for G , M20, C, and A 

G, M20, and C are considered reliable for galaxies that have 
(S/N) per pixel ^ 2.5, elliptical Petrosian radii r e il 0.3", 
and contiguous segmentation maps (as determined by J. 
Lotz's code; Lotz et al. 2004, 2008a). Reliability criteria for 
measurements of A differ only in that the (S/N) per pixel 
must exceed 4.0, instead of 2.5, a lower limit similar to that 
suggested by Lotz et al. (2004). 

3.1.6 Sersic profiles 

The Sersic profile of a galaxy is given by 



E(r) = E e exp < — k 



l/n 

-I - 1 

r, 



(6) 



where n is the Sersic index of the galaxy; n = 1 represents 
a purely exponential light profile (a disc galaxy), and n = 4 
represents a purely de Vaucouleurs light profile (an elliptical 
galaxy). 

Following Haussler et al. (2007) and Griffith et al. (in 
preparation), we use GALFIT Version 2.1c (Peng et al. 2002) 
and the HST/ ACS V band PSF to measure the Sersic in- 
dices of the original and simulated galaxies in our sam- 
ple. For the initial guesses allowed by GALFIT, we use 
the following values: n = 2.5, b/a = 1-ELLIPTICITY, r e = 
0.162xFLUX_RADIUS 1 ' 87 , and mag=MAG_BEST; ELLIPTICITY, 
FLUX_RADIUS, and MAG_BEST are determined by SExtractoi0 
(Bertin & Arnouts 1996). 

Reliable GALFIT results are selected by requiring rel- 
ative uncertainties of less than 15 per cent of the effective 

3 ftp : //ftp . iap . fr /pub /from_users /bertin / sextractor / 
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radii and the Sersic indices (Ar e ^ 0.15xr e ; An ^ 0.15xn), 
GALFIT magnitudes and SExtractor magnitudes agreeing 
to within 1.0 mag (|?TiGALFIT TTlSExtractor < 1.0), and 
Sersic indices 0.2 < n < 1.5 or 1.5 < n < 8. It is also 
required that GALFIT converges on a value and does not 
crash while analyzing the galaxy (Griffith et al., in prepara- 
tion). 

We exclude from our analyses galaxies to which GAL- 
FIT has assigned a Sersic index n = 1.5 because we find 
that an increase of the fraction of light contributed by a 
simulated AGN (particularly for AGN fractions > 30 per 
cent) corresponds to an increase in the number of galax- 
ies having Sersic indices n — 8 and n — 1.5. The former 
is expected, as the PSF used to simulate the AGNs has a 
Sersic index n = 8, but the latter is unexpected. Visual ex- 
amination of the excluded galaxies reveals unusually promi- 
nent nuclear regions, such that we do not expect to observe 
similar galaxies. Therefore, excluding such galaxies does not 
limit the applicability of our analyses to observed AGN host 
galaxies. 

A plausible explanation is that the PSF component of 
these systems dominates the fitting process of GALFIT, and 
thus the galaxy profile is not properly measured; whether 
the Sersic index is assigned n — 8 or n = 1.5 depends on 
the specific characteristics of the galaxy and the profile fit- 
ting (C. Y. Peng, private communication). Although a full 
analysis of the cause behind this is beyond the scope of this 
paper, the convergence at n = 1.5 disappears when we in- 
clude PSF fitting while running GALFIT on a subset (~ 40) 
of the simulated AGN host galaxies. Thus the convergence 
at n — 1.5 may result in part from not fitting for a point 
source. 

3.1.7 Definitions of morphology classifications 

Based on calibrations using more than 200 low-redshift 
galaxies, Lotz et al. (2004) found that the sensitivity of both 
G and M20 to multiple nuclei makes the G-M20 classifica- 
tion technique particularly well suited for the identification 
of interacting galaxies. They also defined three morphology 
classifications for the non-interacting galaxies. Lotz et al. 
(2008a) refined the interacting and non-interacting classifi- 
cations for high-redshift galaxies and additionally suggested 
classifications based on combining G with A, which is also 
commonly combined with C (Conselice 2003). Table 1 sum- 
marizes the morphology classification definitions adopted for 
this paper; the G-M20 and C-A definitions exactly follow 
Lotz et al. (2008a) and Conselice (2003), respectively, while 
the G-A definitions combine the definitions suggested by the 
two references. 

Morphology classifications based on the Sersic indices 
are also summarized in Table 1, and are based on work by, 
e.g., Blanton et al. (2003), Trujillo et al. (2007), and Bell 
(2008). The 'bulge+disc' classification indicates a Sersic in- 
dex between the disc-dominated and bulge-dominated clas- 
sifications. 

These classifications do not always suggest the same 
morphology for a given galaxy. For example, a single galaxy 
may be classified as interacting by G-M20 and G-A, early- 
type by C-A, and bulge+disc by the Sersic index [e.g., Fig- 
ure[2] Row (a)]. Visual morphology classification is therefore 
still an important check on the morphologies based on the 



Table 1. Morphology classification definitions. 'Merger' is a sub- 
category of 'Interacting' that specifically refers to the galaxy coa- 



lescence stage of 


an interaction. 




Morphology 


Definition 




G-M 20 


Interacting 

E/SO/Sa 

Sb-Ir 


G > -0.14 M 20 + 0.33 

G < -0.14 M 20 + 0.33 & G > 0.14 M 20 - 

G < -0.14 M 20 + 0.33 k. G < 0.14 M 20 - 


H0.80 
t-0.80 


G - A 


Interacting 
Non-interacting 


G > -0.4 A + 0.66 or A > 0.35 
G < -0.4 A + 0.66 & A < 0.35 




C - A 


Merger 

Early-type 

Late-type 


A > 0.35 

A < 0.35 & C> 3.5 
A < 0.35 & C < 3.5 




Sersic Index n 


Disc 

Bulge+disc 
Bulge 


0.2 < n < 1.5 
1.5 < n < 2.5 
2.5 < n < 8.0 





location in G/M20/CAS/ Sersic index space. However, the 
focus of the current study is an investigation of the effect 
of an AGN on the automated morphology measurements, so 
we adopt the morphology classifications determined by the 
automated methods. 

3.2 Simulation of an optically visible AGN 

We use the Tiny Tim software (Krist 1993), which was 
specifically designed to model the PSFs for observations 
taken with various HST cameras, to determine appropriate 
V- and /-band PSFs for the AEGIS HST/ ACS observations. 
For the current work, the PSF is contained within a square 
array of 121 pixels x 121 pixels, designed to complement 
the size of postage stamps created from the original galaxy 
images. After selecting the 622 non-AGN galaxies that con- 
stitute our initial sample, the presence of an AGN is simu- 
lated in each galaxy by adding the HST/ ACS PSF to each 
galaxy's V band image. The same PSF is provided to GAL- 
FIT for its determination of the Sersic profiles. However, we 
do not use GALFIT to explicitly fit for the presence of the 
PSF, rather it is used for image convolution to account for 
blurring (Peng et al. 2002). 

As described by, e.g., Kim et al. (2008), the choice of 
PSF can significantly affect galaxy morphology measure- 
ments. Although it is not the goal of the current paper to 
investigate PSF fitting, we do want to understand how the 
choice of PSF may affect our results. Therefore, we also run 
GALFIT providing a PSF based on HST/ ACS observations 
and find that our results do not significantly depend upon 
the input PSF. 

The HST/ ACS images used for the current study have 
been drizzled, but the PSF has not. Because G and M20 
specifically analyse the light within pixels, we want to know 
how this may affect our morphological analyses. To do this, 
we select a star from one of the drizzled ACS images used for 
the current work that is approximately the same size as the 
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PSF. After running the morphology code on both the star 
and the Tiny Tim PSF, we find good agreement between 
the two sets of measurements; the differences between the 
morphology measurements are \AG\ = |AAf2o| = O.Of. 

Taking the SExtractor-derived value MAG_BEST as an ap- 
proximation of the apparent magnitude of each galaxy, the 
V band fluxes of the initial galaxies are calculated as fol- 
lows: f gal = i -o.4x(MAG_BEST-26.486) ) whefe ^ hag unitg 

of [count s/s]. We then determine the scaling factors neces- 
sary to cause the PSF to contribute a fraction %AGN of the 
total V band flux observed from the resulting galaxy /AGN 
system, where fAGN = fgal X [%agn/(1 - %agn)]. Values of 
%agn range from per cent (the initial galaxy) to 50 per 
cent, and increase by 5 per cent increments (%agn = per 
cent, 5 per cent, 10 per cent, 15 per cent, etc.). 

When a galaxy's natural position of maximum inten- 
sity (as defined by SExtractor's X_PEAK and Y_PEAK; Bertin 
& Arnouts 1996) is offset from the galaxy's geometrical cen- 
tre, increasing the brightness of the geometrical centre may 
result in a galaxy image containing two bright nuclei. Al- 
though this could allow insightful tests of the morphology 
measurements, the current goal is to increase the intensity 
of the natural nucleus. Therefore, when the scaled PSFs are 
added to the galaxy images, the centre of the PSF array is 
positioned to coincide with (X_PEAK, Y_PEAK) from the SEx- 
tractor measurements of the original galaxy images. 

Figure [2] shows three of the galaxies in which an AGN 
has been simulated. The original HST/ ACS V band image 
of each galaxy is shown, followed by V band images of the 
galaxy for which 5, 10, 20, and 30 per cent of the observed V 
band light is contributed by an AGN, with the AGN fraction 
increasing to the right. Rows (a), (c), and (e) show the orig- 
inal and altered galaxy images. For rows (b), (d), and (f), 
we overplot the images from the preceding rows with circles 
indicating the original and current Petrosian radii. For all 
of the images, the light is shown on a log scale, using the 
same minimum, maximum, contrast, and bias (0, 40, 2.02, 
and 0.247, respectively). These three galaxies initially meet 
the reliability criteria (§ I3.1.5|) for all of the morphology 
measurements described above, and the figure caption pro- 
vides the G-M20, G-A, C-A, and Sersic profile morphology 
classifications based on the original galaxies. 

Based on the visibility and distinct outlines of the sim- 
ulated AGNs, the AGN shown in row (a) is visually identi- 
fiable at an AGN fraction of 20 per cent, the AGN shown 
in row (e) is visually identifiable at an AGN fraction of only 
10 per cent, and the AGN shown in row (c) appears visually 
identifiable at an AGN fraction of 20-30 per cent. The pres- 
ence of a bright nucleus does not seem out of place for the 
elliptical galaxy shown in row (c), complicating the visual 
distinction between a star forming bulge and an AGN for 
such galaxies. As indicated by the yellow (inner) circles in 
Figure [2] the Petrosian radii decrease significantly for AGN 
fractions exceeding 20-30 per cent. Thus, the size of the Pet- 
rosian radius correlates with the visibility of the AGN; as the 
former decreases, the latter increases. 



3.3 Results 

In this section, we describe the effects of an AGN on the 
elliptical Petrosian radius and individual (e.g., G, C) and 
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Figure 3. Elliptical Petrosian radius r e n as a function of %agn- 
Symbols represent median radii for the full sample of galaxies 
(top), bright galaxies (Mb < —20; middle) and faint galaxies 
(Mb > —20; bottom). Black error bars indicate the standard 
deviations for each %agn- The dotted line indicates the minimum 
required radius, r e n = 0.3", for the morphology measurements 
used in the current study. 



paired (e.g., G-M20) morphology measurements. In the fig- 
ures, the symbol shapes and colours indicate the initial mor- 
phology classifications of the galaxies represented. For ex- 
ample, when presenting the results for G and M20, filled 
red circles represent systems that are initially classified by 
G-M20 as E/SO/Sa galaxies. The number of objects repre- 
sented in each figure is indicated on the individual plots; it 
is found to rapidly decline for %agn > 20 per cent, most 
often due to a shrinking r e n. 

3.3.1 Individual morphology measurements 

The reliability criteria of A, C, G, and M20 measurements 
require that a galaxy's elliptical Petrosian radius, r e ii, ex- 
ceed 0.3" (§ I3.1.5|l . However, as the light observed from a 
galaxy /AGN system is increasingly dominated by the AGN, 
the measured r e ii decreases toward the size of the AGN. 
Figure [3] shows the median elliptical Petrosian radii r e n 
for the entire galaxy sample (top panel), bright galaxies 
(Mb < —20; middle panel), and faint galaxies (Ms > —20; 
bottom panel). At AGN fractions of 25 per cent or more, 
the median radii are similar for all three samples, while at 
AGN fractions of 20 per cent or less, the median radii of the 
brighter galaxies are slightly larger than the median radii of 
the full and faint samples. 

For AGN fractions greater than 20 per cent, the median 
Petrosian radii of galaxies in the full sample are smaller than 
the required 0.3" (dotted line in Figure |3J). This causes a sig- 
nificant decrease in the number of galaxies for which the 
morphology measurements studied here are reliable, which 
has two serious implications. (1) The galaxies that continue 
to provide reliable morphologies at the higher AGN fractions 
are not representative of the general population of AGN host 
galaxies; these are systems for which the extended regions 
are bright enough to successfully compete with the increas- 
ingly bright AGN light in the nuclear regions. The trends 
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Figure 2. Three examples of non-AGN galaxies in which we simulate an AGN represented by a point spread function. For each set, the 
fraction of B band light contributed by the AGN (0, 5, 10, 20, and 30 per cent) increases to the right. The images are shown on a log 
scale ranging from 0-40, with a contrast of 2.02 and a bias of 0.247. Row (a): Galaxy 13009807: z = 0.335; G-M20 '■ interacting; G-A: 
interacting; G-A: early- type; Sersic profile: bulge+disc. Row (b): Same galaxy as row (a); the radius of the white circles is the Petrosian 
radius measured for the initial galaxy, and the radii of the yellow circles are the Petrosian radii measured for the galaxies with PSFs 
added; the bar is 0.5 " in length. Row (c): Galaxy 13018030: 2 = 0.450; G-M20 : E/SO/Sa; G-A: non-interacting; G-A: early-type; Sersic 
profile: bulge-dominated. Row (d): Same galaxy as row (c), with circles as described for row (b). Row (e): Galaxy 13017884: 2 = 0.385; 
G-M20: Sb-Ir; G-A: non-interacting; C-A: late-type; Sersic profile: disc-dominated. Row (f): Same galaxy as row (e), with circles as 
described for row (b). 
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Figure 4. AA as a function of %agn- The symbols indicate the 
median changes AA = Afi na i — Aj n j t ; a i, black error bars indicate 
the standard deviations for each %agn 5 and gray error bars in- 
dicate the maximum and minimum changes at each %agn- Sym- 
bol shapes and colours indicate whether the original galaxies are 
classified by C-A as a merger (black triangles), an early- type (red 
circles), or a late- type (blue squares). The numbers in parentheses 
indicate the initial number of galaxies assigned to each morphol- 
ogy classification, and the number of galaxies having reliable A 
measurements at each %agn is indicated next to the correspond- 
ing symbol. 



4 
2 

-2 

CO 

13 4 

'S 2 
I 

c 

-2 
4 
2 

-2 



54-1 



u 



Merger ( 13) 



•4 4 4 *-{ A 

1 h ; — I 1 T- 



Early-type (120) 
17 



.1.161... 117 1 ! 



97. 



Late 



type (249) 
237, 



*| 49 | 26| 



3 # 



| 22*11 U4j 148| lB j 4 j 



10 15 20 25 30 35 40 
AGN 



Figure 5. AC as a function of %agn- The symbols indicate 
the median changes AC = Cfi na i — C; n i t i a i, and the error bars 
and symbols are as in Figure [4] The numbers in parentheses in- 
dicate the initial number of galaxies assigned to each morphol- 
ogy classification, and the number of galaxies having reliable C 
measurements at each %agn is indicated next to the correspond- 
ing symbol. The measurements either increase as expected or are 
minimally affected for %agn < 20 per cent; above this the mea- 
surements become unreliable. 



shown in the next several figures should be considered with 
this in mind. (2) Figure [3] suggests that most of the AGNs 
in systems for which we do have reliable morphology mea- 
surements contribute little more than 20 per cent of the 
observed optical light. This conclusion is supported by re- 
sults from Pierce et al. (2007), who studied two samples of 
AEGIS AGN host galaxies at redshifts 0.2 < z < 1.2 and 
measured reliable morphologies (G and M20 ) for at least 
80 per cent of their AGN samples. This suggests that the 
majority of AGNs at z ~ 1 contribute less than 20 per cent 
of the optical light observed from the galaxy/ AGN systems. 

Unless the simulated AGN is positioned at a loca- 
tion significantly different from the geometrical centre of 
a galaxy, the change in A as %agn increases is minimal, 
as shown by Figure [H a plot of the median changes AA = 
^4finai — Ai n i t i a i as a function of the AGN fraction. Galaxies 
initially classified as mergers do appear increasingly asym- 
metric as %agn increases to 25 per cent; the decrease for 
higher AGN fractions should indicate that the morphology 
code is focusing on the PSF, instead of the full galaxy 

The focus of the morphology code on the PSF is also 
evident by the decreasing number of galaxy/ AGN systems 
that provide reliable morphology measurements. At an AGN 
fraction of 10 per cent, 31 per cent (191/622) of the systems 
do not have reliable asymmetry measurements; 14 per cent 
(27/191) of these drop-outs are measured to have ellipti- 
cal Petrosian radii r c n < 0.3". For an AGN fraction of 20 
per cent, we find a slightly higher overall drop-out rate (33 
per cent; 207/622), but a significant fraction (72 per cent; 
150/207) of these are are excluded because r c n < 0.3". This 
pattern continues as the AGN fraction increases. 

Figure [5] shows that increasing the AGN fraction gener- 



ally causes an increase or a minimal change in C for %agn < 
20 per cent. As defined in § 13.1.21 C is a ratio of the radii 
enclosing 20 per cent and 80 per cent of a galaxy's light, 
thus concentration measurements should only be reliable for 
%agn < 20 per cent, otherwise the inner radius simply en- 
closes the PSF that has been added to the initial system, and 
the outer radius approaches the same location as the %agn 
increases. For certain host galaxies, the sensitivity of C to 
the defined galaxy centre may provide an alternate explana- 
tion for negative values of AC. This would be particularly 
relevant for AGNs offset from the galactic centre. 

As for the A measurements, the number of reliable 
meausurements of C depends on the fraction of light con- 
tributed by the AGN. Twenty per cent (124/622) of the C 
measurements are deemed unreliable at an AGN fraction of 
10 per cent; we measure elliptical Petrosian radii r e ii < 0.3" 
for 22 per cent (27/124) of these. For an AGN fraction of 
20 per cent, 56 per cent (348/622) of the systems do not 
provide reliable morphology measurements, and 43 per cent 
(150/348) of these are due to small elliptical Petrosian radii. 

Figure [6] shows the median changes AG = Gfl na i — 
Gi n itiai for each of the initial morphology classifications. It 
is clear from this figure that Sb-Ir galaxies are the most sig- 
nificantly affected by the presence of a bright AGN, and 
E/SO/Sa galaxies are the least affected. Several galaxies, 
particularly those initially classified as interacting, exhibit 
a decrease in G as the AGN fraction is increased, instead of 
the expected increase. The disturbed morphologies of some 
of these galaxies, combined with the added PSF, may re- 
sult in an inconsistently defined boundary to the galaxy, 
thereby confusing measurements of G. However, G may be 
more reliable than C for high %agn, in part because G is 
less sensitive to a defined galaxy centre. 

The number of systems for which G is reliable also de- 
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Figure 6. AG as a function of %agn- The symbols indicate the 
median changes AG = G nna i — Gj n itial> and the error bars are as 
in Figure[4] Symbol shapes and colours indicate whether the orig- 
inal galaxies are classified by G-M20 as interacting (black trian- 
gles), E/SO/Sa (red circles), or Sb-Ir (blue squares). The numbers 
in parentheses indicate the initial number of galaxies assigned to 
each morphology classification, and the number of galaxies having 
reliable G measurements at each %agn is indicated next to the 
corresponding symbol. Sb-Ir galaxies are the most significantly 
affected by the presence of a bright AGN, and E/SO/Sa galaxies 
are the least affected; interacting galaxies provide an unexpected 
result with several examples of decreasing G, instead of the ex- 
pected increase. 



pends on the AGN fraction. At an AGN fraction of 10 per 
cent, 22 per cent of the unreliable measurements are caused 
by r cii < 0.3", 22 per cent have low (S/N) per pixel, and 
80 per cent (99/124) of the unreliable systems are deemed 
non-contiguous (§ 13.1.50 ; some systems fail multiple criteria, 
so the sum of the fractions exceeds 100 per cent. In con- 
trast, at an AGN fraction of 20 per cent, elliptical Petrosian 
radii r e u < 0.3" account for 88 per cent of the unreliable 
systems, low (S/N) per pixel account for 15 per cent, and 
non-contiguity accounts for 19 per cent (again, some sys- 
tems fail multiple criteria). Thus, non-contiguous systems 
are most problematic at lower AGN fractions, while small el- 
liptical Petrosian radii are most problematic at higher AGN 
fractions. 

As %agn increases toward 20 per cent, M20 tends to 
decrease for Sb-Ir galaxies and interacting galaxies, as ex- 
pected, while the measured values of M20 for E/SO/Sa galax- 
ies tend to remain unchanged. These effects are demon- 
strated in Figure [7] a plot of the median changes AM20 = 
M20, final — M20, initial- The reliability criteria for M20 are iden- 
tical to those for G, so the explanation for the decrease in the 
number of reliable measurements is the same as described 
in §G£1 

The unexpected increases in the M20 measurements of 
E/SO/Sa galaxies and the large dispersions in AM20 for Sb- 
Ir and interacting galaxies at %agn ~ 20 per cent can be 
at least partially understood by a review of the definition of 
this measurement (£j |3.1.4|) . Of particular importance are the 
locations of the brightest 20 per cent of a galaxy's pixels. If 
an AGN contributes at least 20 per cent of the light from the 



Figure 7. AM20 as a function of %agn- The symbols indicate 
the median changes AA^20 — -^20, final — M 20 ,initiai; error bars arc 
as in Figure|4]and symbol shapes are as in Figure[6] The numbers 
in parentheses indicate the initial number of galaxies assigned to 
each morphology classification, and the number of galaxies having 
reliable M20 measurements at each %agn i s indicated next to the 
corresponding symbol. The median changes for E/SO/Sa galaxies 
are not as strong as the changes experienced by galaxies in the 
other two morphology classifications. 



system, many of the brightest pixels will be associated with 
the AGN, and thus restricted to a particular location. In 
general, this will significantly bias the measurement toward 
lower values of M20 ■ On occasion, the galaxy may have a 
few pixels outside the vicinity of the AGN that are bright 
enough to compete with the AGN light, so the classification 
is not inherently doomed for bright AGNs, but it is very 
likely to be misleading. 

The Sersic index n is a measure of a galaxy's light pro- 
file, an indication of whether the galaxy is better represented 
as a disc (n = 1) or a bulge (n — 4). Increasing the bright- 
ness of a few adjacent pixels would tend to increase the 
Sersic index toward a bulge profile. 

As shown in Figure [S] the presence of an AGN con- 
tributing less than 30 per cent of the rest-frame B band 
light affects n approximately as expected for most of the 
galaxies considered. The median increase is weaker among 
disc-dominated systems than among bulge-dominated or 
bulge+disc systems, though the maximum changes are 
greatest among the disc-dominated systems, suggesting that 
the fraction of AGN host galaxies that are classified as 
disc-dominated should be considered a lower limit to the 
true fraction of disc-dominated AGN hosts. Similarly, the 
fraction of AGN host galaxies classified as bulge-dominated 
should be considered an upper limit. 

Simmons & Urry (2008) measured the change An for 
a large sample of real and simulated galaxies with point 
sources added. For input point source magnitudes AI > 22, 
they found that GALFIT is generally able to recover the 
correct n, with an apparently smaller scatter than what is 
presented here. Differences between our results and those 
of Simmons & Urry (2008) may be partially explained by 
the differences in the sample sizes. While the initial sam- 
ple presented here contains 622 galaxies, that presented by 
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Figure 8. An as a function of %agn- Symbols indicate the 
median values of An = nfi na i — n; n i t j a i, and error bars are as 
in Figure U Symbol shapes and colours indicate whether galax- 
ies are initially classified as bulge-dominated (red circles), disc- 
dominated (blue squares), or bulge+disc (black triangles). The 
numbers in parentheses indicate the initial number of galaxies 
assigned to each morphology classification, and the number of 
galaxies having reliable Sersic indices at each %agn is indicated 
next to the corresponding symbol. The median changes experi- 
enced by initially disc-dominated systems are clearly less than the 
changes experienced by initially bulge+disc or bulge-dominated 
systems. 



Simmons & Urry (2008) contains 12,592 galaxies, provid- 
ing more statistically significant results. In addition, they 
combined results from all galaxy types, which may hide the 
differences shown in Figure [8] 



3.3.2 Combined morphology measurements 

Figure [9] shows the results of morphology classifications 
based on a combination of C, which is expected to increase 
as the AGN fraction increases toward %agn = 20 per cent, 
and A, which is expected to experience minimal change, un- 
less the AGN is offset from the galaxy centre. At %agn = 
10 per cent, more than 90 per cent (106/116) of the galaxies 
initially classified as early-type retain their classification. In 
contrast, only 41 per cent (96/237) of the late-type galax- 
ies are still classified as such, and almost half (46 per cent, 
110/237) are classified as early-type. Galaxies initially classi- 
fied as mergers tend to retain that classification at %agn 
1 per cent; the seven mergers with reliable morphologies 
at %agn = 20 per cent are split between mergers (four), 
early- type (one), and late- type (two). 

The C-A classification method has the advantage that 
few of the early-type galaxies are incorrectly classified as 
the AGN fraction increases. However, the early-type classi- 
fication becomes significantly contaminated with late-type 
galaxies, and C measurements are unreliable for AGN frac- 
tions %agn > 20 per cent. For AGN hosts, the early-type 
classification is rather complete, but impure, while the late- 
type classification is mostly pure, but incomplete; the merger 
classification is contaminated with late-type galaxies and in- 



Figure 9. C-A diagrams of simulated AGN host galaxies. The 
values in the lower right corners of each panel indicate the frac- 
tion of rest-frame B band light contributed by the AGNs repre- 
sented in the panel and the number of host galaxies for which 
we have reliable morphologies; only those systems that meet the 
reliability criteria are included. Solid lines separate the galaxies 
by morphology classification, as indicated on the final panel. The 
symbol shapes and colours indicate the initial C-A morphologies, 
following Figure [4] 



0.4 

c 

a 0.7 



_. 1 . . . 1 . . 


1 ,_l 

A A - 
A 

A 

A - 
A 

A A 

A 

384 J 

0% ! 

; i i i i i i i i i i i i i i. 




i ... i ... i ... i ._ 

-ii ' \ " 

j.yr j 

367 J 

10% : 


^ ' '"' 1 ' V 

j§||s£\*.. 

'■_ "a 


. j a \ 

305 j 

20% ; 


-' 1 '.*t 1 1 1 1 

. 

: Non- 

i interacting 


1 1 1 1 1 1 1 1 1 1^1 1 1 l_ 

Interacting " 

99 1 
30% : 



0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
Asymmetry 



Figure 10. G-A diagrams of simulated AGN host galaxies. The 
layout of this figure follows Figure [9] Galaxies initially classified 
as interacting (non-interacting) by G-A are represented by black 
triangles (red circles). 



complete due to a decrease in the asymmetry of a few galax- 
ies. 

In Figure[TD]we demonstrate the results of pairing G and 
A to classify galaxies as interacting or non-interacting. As 
previously discussed, A changes little with increasing %agn 
unless the AGN is located away from the geometric centre 
of the galaxy, so the main result shown in this figure is due 
to changes in G, which tends to increase with %agn (cf. 
Figure . 

Among galaxies initially classified as interacting by the 
G-A method (black triangles in Figure I10[) . very few are 
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Figure 11. G-M20 diagrams of simulated AGN host galaxies. 
The layout of this figure follows Figure [9] Symbols indicate the 
initial morphologies, following Figure [6] 



Figure 12. G- Sersic index diagrams of simulated AGN host 
galaxies. The layout of this figure follows Figure [9] The symbol 
shapes and colours indicate the initial Sersic profiles, following 
Figure [8] 



misclassified as non-interacting as the fraction of light com- 
ing from the point source increases; as a result, the non- 
interacting classification remains pure. However, as the AGN 
contribution increases, the region indicating interaction suf- 
fers from contamination by misclassified non-interacting 
galaxies. For %agn = 10 per cent, 21 per cent (72/343) of 
the non-interacting galaxies are instead classified as inter- 
acting; for %agn = 20 per cent, this fraction increases to 71 
per cent (204/288). Thus, the fraction of AGN hosts classi- 
fied by G-A as interacting represents an upper limit, and the 
fraction of AGN host galaxies classified as non-interacting 
should be considered a lower limit. 

As a galaxy's nuclear flux increases, the measurement of 
the distribution of light, G, is expected to monotonically in- 
crease because more of the light is distributed among fewer 
of the pixels associated with the galaxy. In contrast, the 
value of M20 is expected to decrease as the fraction of light 
contributed by an active nucleus increases. Increasing the 
brightness of a few adjacent pixels, such as within the nu- 
clear region of a galaxy, tends to decrease the spatial dis- 
tribution of the brightest 20 per cent of the pixels with- 
out changing the overall spatial distribution of the pixels, 
thereby decreasing M20 relative to its initial value. Together, 
these effects are expected to shift G-M20 measurements to- 
ward the E/SO/Sa classification (§ I3.1.7[) . Results from the 
pairing of G with M20 are presented in Figure [TT] 

Sb-Ir galaxies are significantly affected by the presence 
of an AGN. At %agn = 10 per cent, only 46 per cent 
(132/288) of the initially Sb-Ir galaxies that still meet the 
reliability criteria are again classified as Sb-Ir, 46 per cent 
are classified as E/SO/Sa, and 8 per cent are classified as 
interacting. In contrast, 81 per cent (78/96) of the galaxies 
initially classified as E/SO/Sa, and considered to have reli- 
able morphologies, are 'correctly' classified at %agn = 10 
per cent, while 4 per cent (4/96) are classified as Sb-Ir and 
14 per cent (13/96) are classified as interacting. However, at 
%agn = 20 per cent, only 34 per cent (23/68) of the galax- 
ies initially classified as E/SO/Sa are still so classified; due 
to increases in G, most are instead classified as interacting. 



The galaxies initially classified as interacting tend to retain 
that classification. 

Both G and the Sersic index n are expected to increase 
with an increase in the AGN fraction; Figures [12] and [13] de- 
pict the pairing of these two measurements. Symbols in the 
former indicate the original Sersic profiles, while symbols in 
the latter indicate the original G-M20 morphology classifica- 
tions. The number of objects represented in the two figures 
differs because the set of galaxies that initially meet the 
Sersic index reliability criteria does not completely overlap 
with the set of galaxies that initially meet the G-M20 relia- 
bility criteria. 

The addition of a point source tends to increase the 
value of G so that although the galaxy sample presented in 
Figure [12] initially has a large range of G values, the range 
for %agn ^ 20 per cent is rather narrow. The Sersic index 
also tends to increase as the fraction of light coming from 
an added point source is increased, but this trend does not 
apply equally to all galaxies, as indicated by the wide spread 
in Sersic indices at high AGN fractions. This suggests that 
G is more strongly affected than the Sersic index by the 
presence of a point source in a galaxy. 

At %agn = 10 per cent, 58 per cent (230/395) of the 
initially disc-dominated galaxies retain that classification, 
and at %agn = 20 per cent, 51 per cent (133/260) are still 
correctly classified. AGNs also significantly affect the Sersic 
profiles of systems initially classified as Sb-Ir galaxies, as 
demonstrated by Figure [13] As the AGN fraction increases, 
about half of the systems initially classified as Sb-Ir galaxies 
transition from disc-dominated or bulge+disc classifications 
into bulge-dominated classifications. The bulge+disc galaxy 
Sersic indices are more easily increased, with only 21 per 
cent (12/56) still classified as bulge+disc at %agn = 10 per 
cent. As a result of these two trends, the bulge-dominated 
classification is quickly contaminated. The fraction of AGN 
host galaxies classified as bulge-dominated should be consid- 
ered an upper limit, while the reverse is true for the fractions 
of AGN hosts classified as disc-dominated or bulge+disc. 
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Figure 13. G- Sersic index diagrams of simulated AGN host 
galaxies. The layout of this figure follows Figure [9] The sym- 
bol shapes and colours indicate the initial G-M20 morphologies, 
following Figure IS] 



3.4 Summary of morphology results 

With this first experiment, we have investigated the effect 
of an AGN on the morphological measurements G, M20, C, 
A, and the Sersic index. First, we find that the number 
of galaxies deemed to have reliable morphology measure- 
ments decreases significantly for %agn 20 per cent, due 
almost exclusively to an elliptical Petrosian radius smaller 
than r e n = 0.3". 

C and M20 are adversely affected for AGN fractions 
%agn 20 per cent on account of their dependence on 
the central and/or brightest 20 per cent of a galaxy's light. 
In such cases, the morphology measurements are effectively 
measuring the light from the AGN, with minimal input from 
the rest of the galaxy. Thus, C and M20 are expected to 
be unreliable when classifying the morphologies of galaxies 
hosting bright, unobscured AGNs. 

The Gini coefficient of Sb-Ir galaxies is significantly af- 
fected even for an AGN fraction of 10 per cent, with an 
estimated 54 per cent misclassified as E/SO/Sa or interact- 
ing galaxies. For the same AGN fraction, the Sersic index 
is slightly less affected by the presence of an AGN; 42 per 
cent of the disc-dominated host galaxies are misclassified as 
bulgc+disc or bulge-dominated. 

However, not all of the morphology measurements are 
so adversely affected. A experiences little change as the AGN 
fraction increases, unless the location of the AGN is offset 
from the geometrical centre of the galaxy. Thus, A may still 
reliably distinguish between interacting/merging and non- 
inter act ing/ non- merging galaxies . 

The end result is that an AGN contributing more than 
10-20 per cent of a system's light in a given band (here, the 
rest-frame B band) significantly affects most morphology 
measurements. However, disc galaxies hosting AGNs con- 
tributing more than 20 per cent of the B band light may 
be visually identifiable (§ 13. 2[) . alerting researchers to take 
special care with such systems. At the highest AGN frac- 
tions, AGNs hosted by elliptical galaxies appear to have 
sharper edges than typical stellar bulges, thus they may be 



visually identifiable in certain systems. Visual inspection of 
a sub-sample of the simulated galaxies described here sug- 
gests that about half of the AGNs are visually identifiable 
for AGN fractions greater than 20-30 per cent, but for a 
few unusual cases, the AGNs are not clearly distinguishable 
from a nuclear starburst even at the highest AGN fractions. 



4 EXPERIMENT #2: COMBINING AGN AND 
NON-AGN GALAXY SPECTRAL 
TEMPLATES 

For the second experiment, an investigation of the effect of 
an AGN on measured optical and UV-optical colours, a spec- 
tral template representing a quasar is added to spectral tem- 
plates representing three non- AGN galaxies (an elliptical, an 
Sb galaxy, and a starburst galaxy). The quasar template is 
scaled to contribute specified fractions of the rest-frame B 
band flux of the resulting galaxy/AGN systems, consistent 
with the method used for the first experiment (§[3]). We then 
measure the rest-frame optical and UV-optical colours of the 
resulting systems. The goal of this experiment is to deter- 
mine the fraction of the total flux contributed by an AGN 
that is necessary to move a red sequence or green valley 
galaxy into the green valley or the blue cloud, respectively. 



4.1 Spectroscopic simulation of an AGN 

Before combining the spectral templates described in §§ 12.21 
and !2.3l the rest-frame B band flux from the quasar template 
is scaled to contribute specified fractions (5 per cent, 10 
per cent, 15 per cent, ... 50 per cent) of the rest-frame B 
band flux from the resulting galaxy+quasar spectra. After 
adding the scaled quasar spectrum to the galaxy spectra, 
the NUV— r and U — B colours of each of the original and 
new systems is measured, and the fraction of light from the 
AGN spectra that moved a red sequence or green valley 
galaxy into the green valley or blue cloud is determined. 

SynpholQ, a package available for the Image Reduction 
and Analysis Facilitjjf] includes a task named calcphot that 
measures the flux /„ [erg s _1 cm -2 Hz -1 ] in a specified 
passband for an input spectrum. Using this task we measure 
/„ for the GALEX NUV, SDSS r, and Johnson U and B 
passbands from each of the spectral templates. Scaling of 
the AGN templates with respect to the non-AGN galaxy 
templates is based on the B passband fluxes measured from 
each spectral template as follows: 



% 



0AGN 



(7) 



where %agn represents the specified fraction of the total 
B band flux contributed by the AGN spectral template; 
/ s tot/gai/AGN represents the B band flux contributed by the 
combined system, the galaxy, or the AGN; and fx iot , gal , AGN 
represents the NUV, r, or U flux contributed by the com- 
bined system, the galaxy, or the AGN. The NUV— r and 
U — B colours are then calculated as follows: 



4 http:/ /www. stsci.edu/resources/software_hardware/stsdas 
/synphot 



http: / /iraf.noao.edu/iraf-homepage.html 
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Figure 14. U — B vs. NUV— R of the initial galaxy sample de- 
scribed in § 12,11 [The colours represented in this figure come from 
different data (§ I2.1|l than the colours presented in the next few 
figures 12.21 and 12.3b . so the UV-optical colours shown here 
(NUV— R) differ slightly from the UV-optical colours shown in 
the later figures (NUV— r).] Gray dots represent the available 
rest-frame optical and UV-optical colours. Galaxies for which 
UV-optical colours are available but not considered reliable are 
marked with a black square, and those for which UV-optical 
colours are unavailable are represented by black triangles and set 
to have NUV— R = 6.5. Solid lines indicate the divisions between 
blue, green, and red UV-optical colours, as labeled; dashed lines 
indicate approximate divisions between blue (U — B < 0.9) and 
red (U — B > 1.1) optical colours. Compared to the UV-optical 
colours, the optical colours saturate at U — B ~ 1.2. 



NUV - r = -2.5 x log(/ NUV //r); (8) 
U-B = -2.5xlag(f u /f B ). (9) 

This method of calculating the colour is equivalent to 
first converting the fluxes to AB magnitudes [e.g., NUV 
= 2.5 xlog(/Nuv) +48.60] and then calculating the difference 
between the NUV and r or U and B measurements. 

Roughly following Willmer et al. (2006), Nandra et al. 
(2007), and Weiner et al. (2009), we define red and blue op- 
tical colours as U — B > 1.1 and U — B < 0.9, respectively. 
The red and blue regions defined in the references are based 
on a combination of U — B and Mb, but since we do not 
use Mb for the current study, we approximate the divisions 
between regions using only the U — B colour. Blue, green, 
and red UV-optical colours are defined as NUV— r < 3, 
3 <NUV-r < 4.5, and NUV-r > 4.5, respectively, follow- 
ing, e.g., Wyder et al. (2007). The UV-optical green valley 
encompasses the significant population exhibiting colours 
between the UV-optical red sequence and blue cloud; a simi- 
larly significant population is not observed between the opti- 
cal red and blue regions. Figure [14] demonstrates the general 
consistency between the optical and UV-optical colours of 
the initial sample described in § 12.11 as well as the divisions 
between colours. The lack of a significant optical green valley 
population is evident from this figure. 



Figure 15. NUV— r as a function of the fraction of B band flux 
contributed by an AGN. Red circles represent the elliptical tem- 
plate, magenta stars represent the Sb template, and blue squares 
represent the starburst template. Symbols outlined in black rep- 
resent galaxies at %agn = per cent, symbols filled with white 
represent galaxies at %agn = 20 per cent, and symbols filled with 
black represent galaxies at %agn = 40 per cent. Horizontal lines 
indicate the separations between red, green, and blue UV-optical 
colours. The AGN template has a colour NUV-r = -1.09. 



4.2 Results 

Figure [15] shows the change in NUV— r colour as a function 
of the fraction of B band flux contributed by the quasar 
template. The colours of the systems resulting from the com- 
bination of the AGN and non-AGN galaxy templates tend 
toward the colour of the AGN, NUV-r = -1.09, as the 
AGN fraction increases. 

The quasar template brings the initially red elliptical 
galaxy down into the green valley and then the blue cloud 
by AGN fractions of 5 per cent and 20 per cent, respectively. 
The Sb galaxy is also drawn from its initial location in the 
green valley down to the blue cloud when only 10 per cent 
of the B band flux is contributed by the AGN. Although 
the colour of the starburst galaxy is also affected by the 
quasar template, it is blue to begin with, so the effect is less 
significant. 

These results suggest that the UV-optical colours of 
AGN host galaxies may be more significantly affected by 
an active nucleus than previous work has indicated (e.g., 
Kauffmann et al. 2007). However, Kauffmann et al. (2007) 
specifically excluded galaxies exhibiting quasar-like spectra 
(e.g., AGNs typically associated with a central point source), 
so this might be a result that depends very strongly on the 
type of AGN observed. In particular, AGNs that do not 
exhibit a rising blue continuum, such as LINERs and Type- 
2 Seyfert galaxies, are not expected to affect the measured 
host galaxy colours as significantly as the results shown here. 

The AGN template similarly affects the U — B colours 
of the three non-AGN galaxies, as Figure [Tfjl demonstrates. 
Again, we find that the quasar template significantly affects 
measurements of both the elliptical and the Sb galaxy tem- 
plate colours, bringing them down into the blue cloud by 
AGN fractions of 25 per cent 15 per cent, respectively. 

Comparing Figure [16] to results from Nandra et al. 
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Figure 16. U — B as a function of the fraction of B band flux 
contributed by an AGN. Symbols are as for Figure [T5l Horizon- 
tal lines indicate approximate separations between red and blue 
optical colours. The AGN template has a colour U — B = —0.01. 



Figure 17. U — B vs. NUV— R of the galaxy samples and tem- 
plates described in §[2] Small (black and gray) symbols follow Fig- 
ure !14l Large (coloured) symbols follow Figure lTol For most AGN 
fractions, the observed colours appear consistent with colours of 
galaxies not hosting AGNs. 



(2007; see their Figure 1), we conclude that observed AGN 
host galaxies that exhibit colours U — B < 0.5 may be con- 
sistent with optically red galaxies, the light from which has 
been contaminated by AGNs contributing more than 50 per 
cent of the system's rest-frame B band light. Their estimate 
of AGN-caused colour contamination (solid line on their Fig- 
ure 1) is consistent with our result for an AGN fraction of 
25-30 per cent. 

It is clear that the light from AGN host galaxies ex- 
hibiting colours U — B < 0.5 (e.g., Nandra et al. 2007), is 
most likely contaminated by the quasar light. However, it 
is less obvious what to conclude about the AGN hosts lo- 
cated in the upper region of the blue cloud. Pierce et al. 
(2010) address this further by considering the relationships 
between host galaxy aperture colours and the AGN X-ray 
luminosities and hardness ratios. 

Finally, Figure [T7] again shows the optical and UV- 
optical colours of the galaxy sample shown in Figure [T3] 
(small gray and black symbols), but it also includes the AGN 
colour tracks created by combining Figures [TS] and 1161 Re- 
gardless of AGN fraction, the colours of the Sb+AGN tem- 
plate (magenta stars) are similar to the rest-frame colours 
of observed galaxies, while the starburst+AGN template 
colours (blue squares) are consistent with the colours of ob- 
served galaxies for all but the highest AGN fractions. In con- 
trast, the elliptical+AGN template (red circles) is optically 
redder than most of the observed galaxies suggesting that 
such a colour-colour diagram of AGN host galaxies (partic- 
ularly elliptical galaxies) may facilitate the identification of 
colour contamination. Contaminated AGN host galaxies in 
the UV-optical green valley may exhibit unusually red opti- 
cal colours. 



5 COLOURS VS. SERSIC INDICES OF 
SIMULATED AGN HOST GALAXIES 

Observations indicate a general correlation between the 
colour and the morphology of an undisturbed galaxy (e.g., 



Bohm & Wisotzki 2007; Lotz et al. 2008a; Silverman et al. 
2008) , which suggests that unusual colour-morphology com- 
binations among AGN host galaxies (e.g., blue ellipticals) 
may indicate examples of AGN contamination. In addition, 
Lotz et al. (2008b) showed that it is possible to use the recent 
star formation history and morphology of a merging galaxy 
to determine the observed merger stage. This information, 
combined with observations such as X-ray luminosity and 
obscuration, facilitates the study of AGN lifetimes and the 
possible connections between galaxy interactions and the ini- 
tiation of supermassive black hole growth. 

In this section, we present colour-morphology diagrams 
created by matching results from the optical and UV- 
optical colour simulations (§0 to results from the morphol- 
ogy simulations (§ I3.3|l . For each AGN fraction, we match 
colour results from the quasar+elliptical, quasar+Sb, and 
quasar+starburst templates to the median morphologies 
of galaxies initially classified as having bulge-dominated, 
bulge+disc, and disc-dominated Sersic profiles, respectively. 

Figure [18] a plot of Sersic index vs. UV-optical colour, 
indicates a general anti-correlation between colour and 
Sersic index. As the fraction of light contributed by an AGN 
increases to 25 per cent, the Sersic index generally increases, 
while the colour becomes bluer. For AGN fractions exceed- 
ing 25-30 per cent, the Sersic indices of bulge-dominated and 
bulge+disc systems are found to decrease as the colour be- 
comes bluer, suggesting a reliability limit for systems with 
strong point sources. This figure also suggests that AGN 
hosts with Sersic indices n > 2.5 and NUV— r < 2.5 are 
strong candidates for high levels of AGN contamination 
(e.g., %agn > 20 per cent). 

Few X-ray-selected AGNs are observed with colour- 
morphology characteristics similar to those of the strongly 
contaminated systems (e.g., Pierce 2009), suggesting either 
that most X-ray-selected AGNs contribute < 20 per cent 
of the rest-frame B band light detected from the system or 
that the contributed AGN light is significantly reddened, in 
contrast to the quasar template used here. 
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Figure 18. NUV-r vs. Sersic index of non-AGN galaxies with 
AGNs added. Symbols and lines follow Figure 1151 except that 
the symbols also represent galaxies initially classified as bulge- 
dominated (red circles), bulge+disc (magenta stars), or disc- 
dominated (blue squares). Horizontal lines divide the diagram 
into red, green, and blue UV-optical colours, as defined in §|4] 
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Figure 19. U — B vs. Sersic index of non-AGN galaxies with 
AGNs added. Lines and symbols follow Figurc HSl except that hor- 
izontal lines divide the diagram into red and blue optical colours 
(e.g., Willmer et al. 2006). 



At AGN fractions less than 20 per cent, there is reason- 
able concern that some AGN host galaxies are incorrectly 
classified as green valley galaxies and not easily identified 
as being problematic because of the low AGN fraction; thus 
extra caution should be taken when studying such systems. 
Using aperture photometry to separate the inner and outer 
colours of the AGN hosts can help with this, as discussed 
by Ammons et al. (in preparation) and Pierce et al. (2010). 

We also compare the Sersic indices to the U — B colours 
of our simulated systems, as presented in Figure [191 U — B 
appears less sensitive than NUV— r to the contribution from 
an unobscured AGN, in particular for the initially bulge- 
dominated system (red circles). This system falls in the mid- 
dle of the optical 'green valley' at an AGN fraction of 20 
per cent, in contrast to its UV-optical blue cloud classifica- 
tion at the same AGN fraction. However, with respect to 
the green valley, optical colours are less discriminating than 
UV-optical colours, as noted in § 14.21 



6 DISCUSSION AND SUMMARY 

It is reasonable to expect that galaxies hosting unobscured 
AGNs (i.e., those visible as point sources) may experience 
significant contamination of their morphological and colour 
measurements. This paper describes our method for and re- 
sults from testing the reliability of common measures of 
galaxy morphology and colour when an AGN contributes 
a specified fraction of the observed B band light. We em- 
phasize that the work presented here was specifically de- 
signed to investigate the limitations of various morphology 
and colour measurements in the presence of an AGN, in part 
for the sake of better understanding previous results that 
relied on these measurements without explicitly accounting 
for the AGN light. Others have carefully investigated AGN 
host galaxy morphologies and colours, explicitly removing 
the AGN contribution (e.g., Sanchez et al. 2008; Gabor et 



al. 2009), and the reliability of host galaxy/ AGN decompo- 
sitions with respect to PSF variations (Kim et al. 2004). 

For the first experiment (§§ [3] [3T3J , we simulate AGNs, 
represented by a V band PSF, in a sample of AEGIS galaxies 
not known to host AGNs. We first find that increasing the 
fraction of observed light contributed by an AGN decreases 
the measured elliptical Petrosian radius, which in turn de- 
creases the number of galaxies for which we can determine 
reliable morphologies ( § [3j . This decrease is significant for 
AGN fractions exceeding 20 per cent. Two of the morphol- 
ogy measurements, C and M20, are almost certainly unre- 
liable for AGNs that contribute more than 20 per cent of 
the galaxy's B band light. By definition, these two measure- 
ments depend in part on locating either the central or the 
brightest 20 per cent of the flux; if that flux comes directly 
from an AGN, the measurements focus on the AGN and may 
not properly characterize the host galaxy. The presence of an 
AGN also detrimentally affects measurements of the Sersic 
index and G, as expected, though the effects are not inher- 
ently unreliable. In contrast, unless the AGN is offset from 
the geometric centre of the galaxy, A remains generally re- 
liable even for high AGN fractions. As described in § 13.31 
the results from a similar experiment by Simmons & Urry 
(2008) are complementary to the results presented here. 

Gabor et al. (2009) showed that AGN host galaxy con- 
centrations are lower when the AGN contribution is first fit- 
ted with a PSF and then subtracted from the image, consis- 
tent with our results that the addition of a point source tends 
to increase the measured concentration. However, their re- 
sult that the asymmetry of a galaxy may be artificially low 
if the AGN is not properly accounted for does not seem con- 
sistent with the results presented here. Residuals shown by 
Kim et al. (2008) suggest that PSF subtraction may some- 
times increase the apparent asymmetry of the galaxy, un- 
less the PSF is very well matched to the AGN being fitted. 
Though the current work uses only a single PSF for all of 
the simulations, we do not attempt to subtract it from the 
images; in addition, we provide GALFIT with the same PSF 
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that we add to the galaxy images. Thus we do not expect 
to run into problems with this particular issue. 

For the second experiment (§H]|, we add a quasar spec- 
tral template to three non-AGN spectral templates to deter- 
mine how an AGN affects the optical and UV-optical colours 
of its host galaxy. The quasar spectral template is scaled to 
contribute a set of specified fractions of the rest-frame B 
band light emanating from the final galaxy/ AGN systems, 
and the colours of the original and final systems are mea- 
sured and compared. The use of a quasar spectral template 
avoids the complications which would be introduced by the 
use of a spectral template representing an obscured AGN, 
such as the effect of the obscuring gas on the AGN con- 
tribution to the colour of the system. However, the quasar 
template represents an extreme example of an AGN, mean- 
ing that our results regarding the potential colour biases are 
to be considered upper limits. 

Quasar light strongly affects the measured UV-optical 
colours of the simulated quasar host galaxies. A quasar con- 
tributing as little as 5 per cent of the observed B band light 
can move the observed NUV— r colour of an elliptical galaxy 
from the UV-optical red sequence to the green valley, and 
an AGN fraction of 20 per cent moves it into the blue cloud. 
Optical colours are similarly affected. An AGN fraction of 25 
per cent shifts the observed optical colour of the elliptical 
galaxy from the red sequence into the blue cloud. An en- 
couraging caveat is that there are relatively few AGNs that 
have such unobscured quasar spectra, and we expect that 
such systems would be easily identified by visual inspection 
of the galaxy images. 

Sanchez et al. (2004), who studied a sample of 15 opti- 
cally selected, intermediate luminosity, type 1 AGNs, found 
that the colours of the early-type galaxies, which comprised 
much of their sample, tended to be on the blue edge of the 
red sequence, if not in the blue cloud itself. Having first fitted 
the AGNs with PSFs and subtracted them from the images, 
these colours are expected to accurately represent the galax- 
ies without contamination from the AGN light. In contrast, 
the artificially blue colours that result from the combina- 
tion of a QSO spectrum and a quiescent galaxy spectrum, 
as shown here, represent an upper limit to the contamination 
from AGN light. 

By combining the results from our two experiments 
(§©i we nn d that systems for which an AGN contributes 
more than 20 per cent of the B band light appear signifi- 
cantly bluer and more bulge-dominated than the underlying 
galaxy. Therefore, we would expect that AGNs hosted by ap- 
parently blue, bulge-dominated galaxies bias the colour and 
morphology measurements. Pierce (2009) found that such 
systems make up less than half of the AGN population at 
z ~ 1, and many of the AGNs that do fit that description 
have visible point sources. Red, bulge-dominated AGN host 
galaxies are not expected to experience contamination from 
the AGN light. 

Based on an inspection of our simulated galaxies, it ap- 
pears that the point sources in such systems should be visu- 
ally identifiable for high AGN fractions. The results of Pierce 
et al. (2010), in which the authors compare the observed 
nuclear, extended, and integrated optical colours of galax- 
ies hosting X-ray-selected AGNs to a corresponding control 
sample, suggests a low incidence of systems with AGN frac- 
tions %agn ^ 20 per cent because the integrated colours of 



most of the AGN host galaxies correlate much more closely 
with the extended colours than with the nuclear colours. 

Employing a collection of simulations, the current study 
investigates the extent to which various methods used to 
measure galaxy colours and morphologies may be biased by 
the presence of an AGN. This is particularly relevant when 
(1) considering the significance and timing of the influence 
that an AGN may exert on the star formation rate of its 
host galaxy (i.e., are 'green valley' AGN host galaxies intrin- 
sically red?) and (2) using the morphological characteristics 
of AGN host galaxies to investigate relationships between 
interactions and the activation of the central supermassive 
black hole (i.e., could the many AGN host galaxies classi- 
fied as bulge-dominated systems actually have disc-like or 
disturbed morphologies which are hidden because of light 
from the nuclear regions?). 

Our results suggest that the answer to the second ques- 
tion is tentatively 'no', unless there is a visually identifiable 
point source. The answer to the first question is less cer- 
tain, because we find that the colours may be affected by 
even a small AGN fraction. However, the quasar template 
used here represents an extreme (and relatively rare) case of 
an unobscured AGN; an increase in the obscuration of the 
AGN is expected to decrease the colour contamination (cf. 
Figure 1 of Nandra et al. 2007). Therefore, our results are 
not necessarily inconsistent with earlier results (e.g., Kauff- 
mann et al. 2007), which indicate that galaxies hosting the 
more numerous AGNs (e.g., Seyfert galaxies and LINERs) 
are not likely to suffer significant contamination from the 
nuclear regions. 

Measurements of AGN host galaxies that indicate 
bulge-dominated Sersic profiles and blue colours are not nec- 
essarily biased by light from the AGN. Studying the aperture 
colours of AGN host galaxies should help to resolve the ques- 
tion of the origin of the blue colour. Pierce et al. (2010) find 
that X-ray-selected AGNs hosted by UV-optical green valley 
galaxies typically exhibit extended optical colours that are 
similar to the extended colours of a control sample, which 
suggests that most green valley galaxies hosting AGNs are 
intrinsically green and that it is rare for red sequence galax- 
ies to host AGNs contributing more than ~ 10 per cent of 
the B band light. 
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